Effect of Electron-Phonon Interaction Range for a Half-Filled Band in One Dimension 
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We demonstrate that fermion-boson models with nonlocal interactions can be simulated at finite 
band rilling with the continuous-time quantum Monte Carlo method. We apply this method to 
explore the influence of the electron-phonon interaction range for a half- filled band in one dimension, 
covering the full range from the Holstein to the Frohlich regime. The phase diagram contains 
metallic, Peierls, and phase-separated regions. Nonlocal interactions suppress the Peierls instability, 
and thereby lead to almost degenerate power-law exponents for charge and pairing correlations. 
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Introduction. — Electron-phonon interaction has an es- 
sential influence on the properties of many materials pQ. 
It plays a key role for pairing and superconductivity, mass 
renormalization, and charge ordering phenomena. Tak- 
ing into account quantum lattice fluctuations leads to a 
complex, many-body problem. Consequently, theoreti- 
cal studies usually rely on simplified microscopic models. 
A frequently invoked approximation, in particular for nu- 
merical studies, is to consider a completely local electron- 
phonon coupling as in Holstein's molecular-crystal model 
[2]. However, nonlocal interactions are expected to play 
an important role in materials with incomplete screen- 
ing such as quasi-one-dimensional (quasi- ID) organics 
[3]. Long-range interactions, as described by the Frohlich 
model [4 , have been investigated in the context of high- 
temperature superconducting cuprates [5HS], and were 
found to support light polarons and bipolarons [T0HT2] . 

Exact numerical methods have played an important 
role for the understanding of coupled electron-phonon 
systems. Whereas efficient algorithms exist for Holstein- 
type models at arbitrary band filling [T3HT8] , extended 
interactions could so far be addressed only in the empty- 
band limit, see, for example, [TO ] [T91I22] . and [8] for a re- 
view. Consequently, key phenomena such as the Peierls 
instability [23] were out of reach. The latter occurs in 
quasi-lD systems with commensurate fillings, for exam- 
ple in TTF-TCNQ [24], and drives a transition to a 
Peierls insulator with charge-density- wave order. For the 
Holstein model, it is known that quantum lattice fluc- 
tuations can destroy the charge order [25], leading to 
a metal-insulator quantum phase transition at a finite 
value of the electron-phonon coupling strength. 

Unbiased investigations of the effect of nonlocal inter- 
actions at finite band-filling represent a long-standing, 
open problem. In this Letter, we use the continuous- 
time quantum Monte Carlo (CTQMC) method [26] to 
study a model that interpolates between and includes 
the paradigmatic Holstein and Frohlich limits. 

Model. — We consider a Hamiltonian H = Hq + Hi, 
where H = —t Y^(ij)a( c ia c ja + H - c describes ID 



fermions with nearest-neighbor hopping t, and 
Hi = £ (sM^+f 4 2 )-7£ f(r)Qi+r (Ai ~ 1) • (1) 

i i,r 

The first term describes lattice fluctuations in the har- 
monic approximation, with the phonon frequency ujq and 
the stiffness constant K = ujqM. The second term rep- 
resents the electron-phonon interaction, in the form of a 
nonlocal density-displacement coupling, with the density 
operator hi = hi a and h^ = c\ a c i(T . The matrix 
elements are chosen as [9[ [10] [27] 

/fr) = {r 2 + e " r/ ^ 0<r<L/2, (2) 

where the lattice constant a = 1. For £ — ^ oo, this cou- 
pling represents a lattice version of the Frohlich interac- 
tion [27 . More generally, Hi may be viewed as an ex- 
tended Holstein interaction with screening length £; the 
original Holstein model [2 is recovered in the limit £ — >• 0. 
For cjo —> oo, the model maps onto an attractive, general- 
ized Hubbard model. Our method can be applied to any 
coupling which preserves translational invariance. The 
restriction of r is due to periodic boundary conditions. 

Method. — For electron-phonon problems [18 , the 
starting point is the partition function at inverse tem- 
perature P = 1/&bT, 

Z = J V(^^)e- So ^^ J V(q)e~ Sl ^^^ , (3) 

where i/j, i/j are Grassmann fields, and q denotes phonon 
coordinates. The phonons can be integrated out exactly 
[28] , leading to a purely fermionic action with a nonlocal 
(in space and time) interaction 

S[ = -f f drdr'Y, ["ito - mf " 1] • (4) 

ij 

The phonon propagator takes the form Dj J J = F(i — 
j)D(r - t') with F(i - j) = Erf( r + il j)f( r ) and 
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the Holstein propagator D(r — r'). The interaction 
range in space (time) is determined by £ (cjo)- The 
CTQMC method used here is based on an exact expan- 
sion around 7 = [26]. A hybridization expansion algo- 
rithm for electron-phonon impurity problems also exists 
[29] . Monte Carlo updates consist of adding or removing 
single vertices, and flipping auxiliary Ising spins p~8j[26]. 
The numerical effort scales with the cube of the aver- 
age expansion order. Because of the underlying weak- 
coupling expansion, the CTQMC method is particularly 
efficient for problems with small expansion orders, such 
as the Peierls transition in the adiabatic regime ujq jt <C 1 
[30] . Importantly, the method enables us to study the 
many-electron problem [31]. We have verified that it re- 
produces exact results in the Holstein limit £ — >• 0, and 
in the anti-adiabatic limit ujq — >• 00 where the model 
maps onto an extended attractive Hubbard model. 

Results. — We choose a phonon frequency uo^/t = 0.5. 
In this regime, representative of many materials, neither 
static nor instantaneous approximations are valid, and 
numerical simulations are essential. The dimensionless 
ratio A = e p /2t, where e p is the polaron binding energy 
in the atomic limit (t = 0) [27] and 2t is half of the free 
bandwidth, allows us to compare different £ at the same 
effective coupling strength. We use t as the energy unit, 
and set h = M = 1. All results are for a half- filled band. 

The phase diagram as a function of £ and A, obtained 
from CTQMC simulations with up to L = 42 sites, is 
shown in Fig. [I] Because the Holstein model is recov- 
ered for £ = 0, its previously studied metallic and Peierls 
insulating phases p~5j [32] [33] smoothly extend to £ > 0. 
However, we observe a significant ^-dependence for small 
£, and saturation for larger values. Additionally, for suf- 
ficiently large £ and A, we find a region of phase sep- 
aration or charge segregation. The metallic phase and 
the metal-insulator transition extend all the way to the 
Frohlich limit £ = 00 [A^ 1 = 0.48(2)]; see also Fig. [3} For 
the range of A shown in Fig.[l] phase separation is absent 
for £ < 2, and becomes more favorable with increasing £. 

The different phases can be characterized by the den- 
sity correlator S p (r) = ((n r — l)(no — 1)) and the density 
structure factor S p (q). In the metallic phase [Fig. |2ja)], 
S p (r) shows a power-law decay of 2kp correlations (with 
exponent K p ) as expected from bosonization, and is lin- 
ear for q — ^ 0. Together with exponentially suppressed 
spin correlations (not shown), these findings are consis- 
tent with a bipolaronic Luther-Emery phase [34[ [35] . 

The Peierls state exhibits quasi-long-range 2/cp density 
correlations [Fig. |2^a)], corresponding to charge-density- 
wave order at T = with two electrons of opposite spin 
forming bipolarons on every other site. The phase bound- 
ary for the Peierls transition can be determined from the 
staggered charge susceptibility [32] 
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FIG. 1. (Color online) Phase diagram as a function of inter- 
action range £ and electron-phonon coupling strength A. The 
regions correspond to a metal, a Peierls insulator, and phase 
separation. The metallic and Peierls phases (and presumably 
also phase separation) extend to the Frohlich limit £ = 00. 
Lines are guides to the eye. Here uoo/t = 0.5. 



For fixed /3/L, xpM/L is universal at the critical point, 
and the crossing of curves for different L gives A^ 1 ; for 
example, we have A^ 1 = 0.33(1) for £ = 2 in Fig. [2^c). 
As for the Holstein model [32] [36] , the Peierls transition 
is expected to be of the Kosterlitz-Thouless type also for 
£ > 0. Extended interactions (£ > 0) promote metallic 
behavior by dissociating the onsite bipolarons predomi- 
nant in the Holstein regime, and we see a Peierls insulator 
to metal transition as a function of £ [Fig. ^a)]. By the 
same mechanism, the critical coupling for the transition 
to the Peierls insulator, A^ 1 , increases with increasing £, 
see Fig. [I] The onset of charge order is also reflected in 
the divergent q = 2kp peak in S p (q), see Fig. [2jlo) . 

Phase separation as a result of the phonon-induced at- 
traction manifests itself as a peak at small q in S p (q), 
as shown in Fig. [5|b) . In the phase-separated region of 
Fig. [I] the quantity 7rS p (qi)/qi (with q\ — 2ir/L) — whose 
thermodynamic limit is related to K p in a Luttinger 
liquid — diverges with system size. In the Peierls phase 
K p = 0, as verified on very large systems [32]. Formally, 
K p = 00, reflecting phase separation, implies a divergent 
compressibility [34]. The divergence of ^S p {q\)/q\ in the 
phase-separated region is shown for £ = 6 in Fig. |2|d), 
and we deduce a critical value of A^ s = 0.55(2). There 
are two possible scenarios for the transition from the 
Peierls to the phase separated region. A continuous tran- 
sition would imply a melting of the charge-ordered Peierls 
state before the formation of multipolaron droplets, al- 
lowing for an intervening (narrow) metallic region with 
finite K p . Alternatively, the insulator-insulator transi- 
tion could be of first order. Evidence for the latter pos- 
sibility comes from the occurrence of metastable config- 
urations and hysteresis at low temperatures in our sim- 
ulations. The occurrence of phase separation in models 
with long-range electron-phonon coupling had been sug- 
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FIG. 2. (Color online) (a) Density- density correlation function S p (r), for A = 0.325 and fit — L — 22, in the Peierls insulator 
(PI, £ = 0.5) and the metallic phase (M, £ = 5). (b) Structure factor S p (q) at £ = 6 in the three phases of Fig. [I] We used 
L = 22, and fi = L (L/2) for A = 0.2, 0.45 (0.55). (c) Scaling of the charge susceptibility x P (7r) [Eq. |5|] at £ = 2, defining the 
critical point A^ 1 = 0.33(1) of the Peierls transition. Here fi — L. (d) Finite-size scaling of 7rS p (qi)Jqi = LS p (2tt/ L)/2 with 
fi — L/2. The divergence for A > 0.55 indicates phase separation (PS). All results are for uoo/t = 0.5. 



gested before [20] and observed in analytical work [6] [7] ; 
for short-range interactions, it is suppressed by the ab- 
sence of bound triplet states [20] , but may occur in the 
vicinity of a Mott transition [37 j. Phase separation is 
expected to be suppressed in the presence of additional 
long-range electron-electron interaction [38] . 

In combination with analytical continuation [39], we 
can calculate the single-particle spectral function 

A(k,w) = | £ I <*l \i) I V 3 * 8 + - w ) ' 

(6) 

where \i) is an eigenstate with energy Ei, and Aji = 
Ej—Ei. Results in the metallic phase (A = 0.2) are shown 
in Fig. [3] for the extreme Holstein (£ = 0.1) and Frohlich 
limits (£ = 00). The locus of spectral weight follows the 
free band dispersion, —2t cos k. The exponentially small 
Luther-Emery spin gap is not resolved for the parame- 
ters chosen, and the spectrum is particle-hole symmet- 
ric. Excitations are sharp inside the coherent interval 
[— cjoj^o]? whereas they are substantially broadened as a 
result of multiphonon processes at higher energies [40] . 
As expected for our ID model, the spectrum agrees well 
with the exact bosonization result [41], which predicts 
a hybridization of the spin, charge and phonon modes, 



although spin-charge separation is not visible due to the 
weak coupling and small uj$ [42]. Comparing Figs. [3ja) 
and (b) we see that in contrast to previous work on one 
and two electrons [SJ [lOj [19 - 22 , the impact of the inter- 
action range is remarkably small. This important char- 
acteristic of the many-electron case can be related to the 
absence of significant polaron and bipolaron effects in the 
metallic phase of Fig. [I] 

Figures [3Ja),(b) show A(k,u) in the Peierls phase 
(A = 0.4). In the Holstein regime (£ = 0.1), the spec- 
trum consists of two sets of features. The cosine band 
seen in Fig. [3] has acquired a gap at the Fermi level and 
reveals additional, backfolded shadow bands as a result of 
dimerization and the corresponding doubling of the unit 
cell [43]. These signatures are labeled (1) and (1') in 
Fig. [4ja), respectively. In addition, we find lower-energy 
excitations labeled (2) corresponding to bound soliton- 
antisoliton pairs or, equivalently, polarons [30, 44 , which 
are absent in a homogeneous mean-field solution that 
captures only (1) and (1') [30]. The soliton dispersion 
indicates a mass larger than the electron mass. Because 
their energy at the Fermi level is lower than the Peierls 
gap, doping would lead to the formation of solitons [45] . 

Increasing the interaction range from £ = 0.1 to £ = 4 
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FIG. 3. (Color online) Single-particle spectral function 
A(k,cj) [Eq. in the metallic phase at A = 0.2 for (a) 
f = 0.1, (b) f = oo. Here cj /t = 0.5, f3t = L = 30. The 
dashed lines indicate the Fermi level (uj — 0) and uj — zLcjq- 
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FIG. 4. (Color online) (a),(b) Single-particle spectral func- 
tion A(k,uj) in the Peierls phase at A = 0.4 for (a) £ = 0.1, 
(b) £ = 4. (c) Dynamical charge structure factor [Eq. ^] for 
the same parameters as in (b). Here uJo/t = 0.5, f3t = L = 22. 
The dashed lines indicate (a),(b) the Fermi level and (c) ujo. 
The labels (1), (1'), and (2) in (a) are explained in the text. 



drives the system into the vicinity of the metal-insulator 
transition, see Fig. [T] This is reflected by a much smaller 
gap at the Fermi level, reduced spectral weight of the po- 
laron excitations, and the suppression of shadow bands. 
Consequently, the spectral function becomes quite simi- 
lar (but not identical) to that shown in Fig. [3| and illus- 
trates the continuous evolution of A(k,uj) across A^ 1 . 

A hallmark feature of the Peierls state is phonon soft- 
ening at q = which is visible [30, 46 in the dynami- 




FIG. 5. (Color online) Charge and pairing correlations at 
distance r = L/2 for (a) £ = 0.1, (b) £ = oo. Here UJo/t = 0.5 
and A = 0.2. Lines are fits to a power law f(r) = Ar -71 . 

cal charge structure factor [p q = e iqr (h r — (h r ))/y/L] 
S p (q,e>) = | £ I <*| \j) \ 2 e-^5(A Si - w) . (7) 

The results in Fig. gc) for A = 0.4, £ = 4 [the same 
parameters as in Fig. gb)] reveal a clear signature of 
the renormalized phonon dispersion. The spectrum is 
dominated by the soft phonon mode at uj = 0, q = 2/cp. 
Furthermore, we observe a continuum of particle- hole ex- 
citations, and a charge gap at long wavelengths (q — > 0). 

Finally, we consider the interaction-range effect on the 
competition of charge and pairing correlations. Figure [5] 
shows results for S p (r) and the s-wave pair correlator 
P(r) = (A+A ) (with At = 4 t c^) in the metallic 
phase. For £ = 0.1, Fig. [5|a) reflects the dominance of 
charge correlations previously observed for the Holstein 
model [47]. However, Fig. [5jb) reveals that long-range 
electron-phonon interaction (here £ = oo) suppresses 
charge correlations and thereby results in almost identical 
power-law exponents in both channels. Additional short- 
range electron-electron repulsion is expected to further 
suppress onsite bipolaron formation — promoting nonlo- 
cal pairing — and to lead to a more general phase diagram 
with Mott, metallic, Peierls and phase-separated ground 
states. It will be interesting to explore this issue further, 
both at and away from half filling, also in the light of a 
recently reported phase of correlated singlets [48] . 

Conclusions. — We used exact Monte Carlo simulations 
to study many-electron systems with nonlocal and even 
long-range electron-phonon interaction. Compared to 
the Holstein model, extended interactions suppress the 
Peierls instability — making pairing more favorable — and 
can lead to phase separation. Key implications for 
materials modeling are that interactions of finite but 
small range are well described by Holstein- type models, 
whereas long-range interactions can have substantial ef- 
fects on the balance of pairing and charge correlations. 
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